Atmos. Chem. Phys., 13, 9577-9595, 2013 
www.atmos-chem-phys.net/13/9577/2013/ 
doi: 10. 5 194/acp- 13-9577-20 13 
© Author(s) 2013. CC Attribution 3.0 License. 



Atmospheric § 
Chemistry | (Rfcj 
and Physics | 


Heterogeneous formation of polar stratospheric clouds - Part 1: 
Nucleation of nitric acid trihydrate (NAT) 

C. R. Hoyle 12 , 1. Engel 2 *, B. P. Luo 2 , M. C. Pitts 3 , L. R. Poole 4 , J.-U. Groofi 5 , and T. Peter 2 

'Laboratory of Atmospheric Chemistry, Paul Scherrer Institut, Villigen, Switzerland 
institute for Atmospheric and Climate Science, ETH Zurich, Zurich, Switzerland 
3 NASA Langley Research Center, Hampton, Virginia 23681, USA 
4 Science Systems and Applications, Incorporated, Hampton, Virginia 23666, USA 

5 Institut fur Energie- und Klimaforschung - Stratosphare (IEK-7), Forschungszentrum Jiilich, Jiilich, Germany 
now at: Institut fur Energie- und Klimaforschung - Stratosphare (IEK-7), Forschungszentrum Jiilich, Jiilich, Germany 

Correspondence to: C. R. Hoyle (christopher.hoyle@env.ethz.ch) 

Received: 14 March 2013 - Published in Atmos. Chem. Phys. Discuss.: 22 March 2013 
Revised: 16 August 2013 - Accepted: 21 August 2013 - Published: 27 September 2013 


Abstract. Satellite-based observations during the Arctic 
winter of 2009/2010 provide firm evidence that, in contrast 
to the current understanding, the nucleation of nitric acid tri- 
hydrate (NAT) in the polar stratosphere does not only oc- 
cur on preexisting ice particles. In order to explain the NAT 
clouds observed over the Arctic in mid-December 2009, a 
heterogeneous nucleation mechanism is required, occurring 
via immersion freezing on the surface of solid particles, 
likely of meteoritic origin. For the first time, a detailed mi- 
crophysical modelling of this NAT formation pathway has 
been carried out. Heterogeneous NAT formation was cal- 
culated along more than sixty thousand trajectories, end- 
ing at Cloud Aerosol Lidar with Orthogonal Polarization 
(CALIOP) observation points. Comparing the optical prop- 
erties of the modelled NAT with these observations enabled 
a thorough validation of a newly developed NAT nucleation 
parameterisation, which has been built into the Zurich Opti- 
cal and Microphysical box Model (ZOMM). The parameter- 
isation is based on active site theory, is simple to implement 
in models and provides substantial advantages over previous 
approaches which involved a constant rate of NAT nucleation 
in a given volume of air. It is shown that the new method is 
capable of reproducing observed polar stratospheric clouds 
(PSCs) very well, despite the varied conditions experienced 
by air parcels travelling along the different trajectories. In a 
companion paper, ZOMM is applied to a later period of the 
winter, when ice PSCs are also present, and it is shown that 


the observed PSCs are also represented extremely well under 
these conditions. 

1 Introduction 

The discovery of the Antarctic ozone hole (Farman et al., 
1985) was rapidly followed by the recognition that polar 
stratospheric clouds (PSCs), by providing the necessary sur- 
face for heterogeneous halogen activation reactions, played 
a central role in the observed ozone depletion (e.g. Solomon 
et ah, 1986). Solid PSC particles were detected at temper- 
atures well above the frost point, and it was proposed that 
these particles contained HNO 3 (Toon et ah, 1986; Crutzen 
and Arnold, 1986). Later, experimental work suggested that 
the composition could be that of nitric acid trihydrate (NAT) 
(Hanson and Mauersberger, 1988), and, more recently, obser- 
vational work, such as the balloon-borne mass spectrometry 
measurements reported on by Voigt et ah (2000), have shown 
that solid PSC particles observed above the frost point do in- 
deed have H 2 O and HNO 3 stoichiometries consistent with a 
NAT composition. 

The formation of NAT particles sequesters reactive nitro- 
gen compounds from the gas phase, and because these par- 
ticles can grow to the large sizes necessary for rapid sedi- 
mentation (e.g. Salawitch et ah, 1988), they cause denitrifi- 
cation of the polar stratosphere. While it is mostly the liq- 
uid particles in the polar stratosphere which are responsible 
for the heterogeneous reactions leading to the activation of 
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the halogens (e.g. Molina et al., 1993; Carslaw et al., 1994; 
Ravishankara and Hanson, 1996), the observed dimensions 
of polar ozone loss are only reached because of this denitri- 
fication, which prevents deactivation of the halogens. 

Despite the importance of NAT particles in determining 
polar ozone loss, an exact formation mechanism of these 
particles has remained elusive (e.g. Peter and Groofi, 2012). 
Initially it was proposed that NAT could crystallise out ho- 
mogeneously in liquid stratospheric aerosol droplets (Molina 
et al., 1993; Beyer et al., 1994; Iraci et al., 1994). Later, 
however, laboratory experiments, using 1 cm 3 solution vol- 
umes, showed that the homogeneous nucleation rates of 
NAT in liquid ternary (H2SO4/HNO3/H2O) solution droplets 
were far too low to lead to the observed formation of such 
PSCs in the polar stratosphere (Koop et al., 1995). The up- 
per bound that Koop et al. (1995) determined for homoge- 
neous NAT nucleation corresponds to an air-volume-based 
rate of /^“^’ vol,uppel = 1.8 x 10” 10 cm” 3 airh” 1 , assuming 
typical stratospheric aerosol particle sizes and number den- 
sities. This results in NAT particle number densities lower 
than 10” 7 cm -3 when an air mass is cooled below 7 nat, the 
temperature below which NAT is thermodynamically stable, 
for one week, a period conceivable for the Arctic. Number 
densities of NAT would still remain below 10” 6 cm -3 , even 
if air was below Tnat for 4 months, as might occur in the 
Antarctic. Such low number densities neither result in effec- 
tive denitrification, nor are they measurable. Subsequent ex- 
periments by the same authors (Knopf et al., 2002), using 
microlitre size solution volumes, determined an upper bound 
of /^™' vol ' upper = 3 x 10” 10 cm -3 airh” 1 , confirming the re- 
sults of Koop et al. (1995). 

Later, Biermann et al. (1996) showed, based on bulk sam- 
ples of ternary solution containing solid meteoritic nuclei, 
that these nuclei could at best have 1/e freezing times of 
around 20 months. While these authors aimed at excluding 
meteoritic material as responsible for dense ( n >0.1 cm” 3 ) 
NAT clouds, their derived upper bound corresponds to an air- 
volume -based rate of 7^ T vol ' upper = 7 x 10” 4 cm” 3 airh” 1 
(assuming a stratospheric smoke surface area density of 
5 pm 2 cm” 3 ), which is certainly high enough to produce ob- 
servable NAT clouds. 

In contrast, the nucleation of NAT on ice is very certain 
and has been demonstrated in the laboratory and in the field. 
Through a combination of observations and microphysical 
modelling, Carslaw et al. (1998) and Wirth et al. (1999) 
showed that NAT can form on ice particles in the atmosphere, 
suggesting that nucleation of NAT on ice is very efficient 
and most likely occurs via deposition nucleation on surfaces 
of ice which are not covered by a liquid coating. Via the 
so-called “mother cloud/NAT-rock mechanism”, Fueglistaler 
et al. (2002) developed a potential mechanism by which NAT 
particles initiated by patchy mountain wave ice clouds are 
distributed vortex-wide. 
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However, a study by Pagan et al. (2004) showed that dur- 
ing early December 1999, while the Arctic stratosphere was 
too warm for synoptic ice formation, there was also no ice 
PSC formation in mountain waves, ruling out nucleation on 
ice as a formation pathway for the NAT clouds observed dur- 
ing this time period. Drdla et al. (2003) presented modelling 
work which also suggested that a NAT formation mechanism 
above the frost point is necessary to explain observed denitri- 
fication in the 1999-2000 Arctic winter. Several other stud- 
ies also exist, where constant NAT nucleation rates have been 
applied to model NAT formation irrespective of the presence 
or absence of ice. Using balloon-based measurements car- 
ried out in the Arctic in December 2002, Larsen et al. (2004) 
found that observed PSC properties can be reproduced by ap- 
plying the surface nucleation rate of Tabazadeh et al. (2002), 
reduced by a factor of 10-20, to calculate NAD (nitric acid 
dihydrate) formation from liquid aerosol, followed by instan- 
taneous transformation of NAD to NAT. The rate they de- 
termined corresponds to an air-volume-based rate of 7/^- 
= 2.5 x 10“ 5 cm” 3 airh” 1 . The air masses which they ob- 
served had not been exposed to synoptic temperatures be- 
low the frost point, and there was no mountain wave activity; 
therefore NAT formation on ice could be ruled out, however 
heterogeneous nucleation on other types of solid particles 
could not be excluded. 

More recently, Voigt et al. (2005) reported aircraft ob- 
servations of low number densities of NAT particles in 
air masses which had been at T > Tnat + 10 K some 24- 
36 h before the measurement and then cooled to 7nat > 
T > Tice + 3 K for a period of less than 20 h before the 
measurement. In order to explain these observations, they 
suggested that NAT may also nucleate heterogeneously on 
the surface of meteoritic dust particles. Because the time 
below Tnat could be well constrained (20 h), they were 
able to estimate a volume-based NAT nucleation rate of 
^NAT™* = 8 x 10” 6 cm” 3 airh” 1 , reproducing their observa- 
tion of NAT number densities of 1.6 x 10” 4 cm” 3 . While 
this rate reproduced the observed number densities, because 
of the limited measurement data, it was not possible to de- 
termine how the NAT nucleation rate increases as temper- 
ature decreases below 7\at- it was argued by Voigt et al. 
(2005) that the heterogeneous NAT nucleation rates on mete- 
oritic material measured by Biermann et al. (1996), while too 
low to be responsible for dense NAT clouds, might actually 
be sufficient to produce low NAT number densities within 
a few hours. Since the nucleation rates published by Bier- 
mann et al. (1996) were upper limits and the real rates might 
be considerably smaller, the suitability of meteoritic dust re- 
mains an open issue, which calls for further investigation. 

Also for the Arctic winter of 2002/2003, Groofi et al. 
(2005) found that using a constant volume -based rate of 
J^j = 8 x 10” 6 cm” 3 airh” 1 in CLaMS (Chemical La- 
grangian Model of the Stratosphere), led to the best agree- 
ment between modelled and observed PSC properties and 
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Table 1 . NAT nucleation rates from previous studies, as described 
in the text. All rates except for /(Snat) are air-volume-based, and 
lead to NAT nucleation as soon as the temperature drops below 
?NAT- 


Nucleation rate 

Reference 

^hom.vol, upper = j g x 10 -10 cm -3 airh -l 

Koop et al. (1995) 

^hom.vol, upper = 3 Q x jq -10 cm -3 ajr h -l 

Knopf et al. (2002) 

^ecvokupper = ? Q x jq -4 cm -3 air h -l 

Biermann et al. (1996) 

■^NAT = 2 '9 x 10 -6 cm~ 3 air h _1 

Carslaw et al. (2002) 

-^NAT = 2 ' 3 x I* 3 ” 5 cm~ 3 air h -1 

Larsen et al. (2004) 

= 8.0 x 10“ 6 cm -3 air h -1 

Voigt et al. (2005) 

/^J T = 8.0 x 10 -6 cm -3 air h _1 

GrooB et al. (2005) 

^NAT™ 1 = 90 X 10_6 cm_3 air h_1 

This work 1 

^NAT = /< S NAT) 

This work 2 


1 The volume-based nucleation rate used in this work, for comparison with the newly 
developed parameterisation. 

- The newly developed parameterisation presented here is not air- volume-based, 
and allows for a saturation dependence of the NAT nucleation rate. 


denitrification throughout the winter. A domain-filling, La- 
grangian particle trajectory model was used by Carslaw et al. 
(2002) to simulate the characteristics of NAT particles in the 
Arctic stratosphere in January-March 2000. They also used 
a simple volume -based NAT nucleation rate, and found that 
the observed particle sizes were well reproduced when a rate 
of J^ T = 2.9 x 10 -6 cm -3 airh -1 was used. While the ap- 
parent patchiness of wave ice clouds was first taken as an 
argument against their effectiveness in NAT nucleation, later 
more detailed vortex-wide simulations by the same authors 
showed that this mechanism could potentially indeed explain 
up to 80 % of the observed denitrification (Mann et al., 2005). 
Davies et al. (2005) simulated Arctic denitrification using a 
combination of a Lagrangian particle model and the chem- 
istry transport model SLIMCAT. They compared the mod- 
elled denitrification with HNO3 measurements from satellite- 
and aircraft-based instruments, and found reasonable agree- 
ment for several winters, when using the same NAT nucle- 
ation rate as Carslaw et al. (2002). However, for the win- 
ter of 1994/1995, they found that in order to reproduce the 
observed denitrification, this nucleation rate needs to be in- 
creased by a factor of 4, highlighting the need for physically 
based nucleation rates. The various air-volume-based NAT 
nucleation rates from previous studies are summarised in Ta- 
ble 1. 

During December 2009, space-borne lidar measurements 
conducted from CAL1PSO (Cloud-Aerosol Lidar and In- 
frared Pathfinder Satellite Observations), showed a situation 
where thin and patchy liquid-NAT mixed clouds of unprece- 
dented dimensions (typically 1000 x 2000 km 2 ) were present 
in air masses in the Arctic, despite these air masses not hav- 


ing experienced temperatures below the frost point (Pitts 
et al., 2011). In the present study, this observational data set 
is utilised to develop a physically based parameterisation for 
NAT nucleation on solid, non-ice particles. The advantage of 
using such satellite -based lidar observations is that they cover 
a large area of the Arctic vortex, and when back trajectories 
are calculated starting from the observations, a range of air 
parcel histories are found, allowing the development of a ro- 
bust parameterisation, which represents NAT formation in air 
reaching different supersaturations with respect to NAT, for 
different lengths of time. 

In the next section, the observational data set is described, 
as well as the models which were applied, and the new NAT 
nucleation parameterisation. In Sect. 3, the results of mod- 
elling NAT formation with the new parameterisation are pre- 
sented and compared with a constant nucleation rate ap- 
proach, and in Sect. 4 we provide a summary and conclu- 
sions of the study. In a companion paper (Part 2; Engel et al., 
2013) the role of heterogeneous nucleation in ice formation 
is discussed. 


2 Observational data and model description 
2.1 The CALIPSO observations 

The CALIPSO satellite was launched in April 2006, and in- 
cludes the Cloud Aerosol Lidar with Orthogonal Polariza- 
tion (CALIOP) as its principal instrument. CALIOP mea- 
sures backscatter at two wavelengths, 1064nm and 532 nm, 
using the 532 nm signal to measure the depolarisation of the 
backscattered radiation relative to the plane-polarised light 
of the outgoing laser beam. With its high-latitude coverage 
up to 82°, CALIOP has provided extensive observations of 
PSCs in both hemispheres since its launch. 

In order to derive and test the new NAT nucleation param- 
eterisation, the CALIOP observations made over the Arctic 
in December 2009 are used. During the second half of De- 
cember 2009, stratospheric temperatures were slightly below 
7nat and patchy areas of PSCs were detected on a daily ba- 
sis. The vortex average temperatures were several degrees 
above 7f ros t (Dornbrack et al., 2012), and, additionally, there 
were no observations of ice clouds during this period. In 
Fig. 1, the minimum temperatures with respect to the frost 
point, along all trajectories used in the analysis here (as de- 
scribed in Sects. 2.3 and 3), are shown. Even the lowest tem- 
peratures along these trajectories are 2.5 K above the frost 
point, while the peak of the distribution is at approximately 
5 K above the frost point. Considering that temperatures ap- 
proximately 3 K below the frost point are required to homo- 
geneously nucleate ice, even accounting for small-scale tem- 
perature fluctuations, the probability of undetected ice clouds 
having occurred is extremely low. This set of observations 
and model calculations therefore supports the existence of a 
NAT nucleation mechanism which does not require an ice 
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T-T frost [K] 


Fig. 1. A distribution showing the minimum temperature reached 
along each trajectory, with respect to the frost point. Trajectories 
were calculated back in time from points along the CALIPSO orbit 
path, as described in Sect. 2.3. All trajectories used to model the 
PSC properties along all orbits in Sect. 3 are included in the figure. 

surface. On 31 December, mountain wave ice clouds were 
first observed, followed by a dramatic increase in observa- 
tions of NAT-containing clouds, marking the end of the pe- 
riod when the NAT PSCs can be considered to have formed 
without the presence of ice particles. 

In the analysis here, we make use of the PSC classification 
system of Pitts et al. (2009, 2011), which assigns the detected 
clouds to one of six categories (STS, Mixl, Mix2, Mix2-enh, 
ice, wave ice), depending on the backscatter and depolari- 
sation signals. Though not identical, this classification sys- 
tem is similar to the previously used PSC types, namely lb 
(STS), la (Mixl, Mix2), la-enh (Mix2-enh), and 2 (ice and 
wave ice). For the NAT classes, Mixl includes NAT num- 
ber densities of less than approximately 10 -3 cm -3 , Mix2 in- 
cludes NAT number densities between about 10 -3 cm -3 and 
0. 1 cm -3 , while Mix2-enh includes NAT number densities of 
above 0.1 cm -3 . 

The PSC detection limit of CALIOP is altitude dependent 
(Pitts et al., 2009), and we treat the modelled data in a similar 
way. If the modelled or observed values at a particular point 
are below the detection limit, then a larger horizontal averag- 
ing distance is applied. Thus, for the observations, if no PSC 
is detected over 5 km, the horizontal resolution is decreased 
(i.e. the averaging distance is increased) to 15km, then 45km 
and 135km. The highest resolution for the modelled data is 
25 km (the resolution of the grid from which the trajectories 
used by the model were started), and if no PSCs are detected 
at a point then the averaging distance is increased to 125 km. 

2.2 The microphysieal box model 

The Zurich Optical and Microphysical box Model (ZOMM) 
used in this study has been applied to simulate PSCs (e.g. 
Meilinger et al., 1995; Luo et al., 2003) as well as cirrus 


clouds (e.g. Hoyle et al., 2005; Brabec et al., 2012). In the 
latter two studies, a slightly different version of the model 
was used; however the only significant difference is that the 
stratospheric version of the model simulates the evolution 
of a population of ternary solution aerosol droplets (H 2 O, 
H 2 SO 4 , HNO 3 ) and their interaction with the gas phase, 
whereas the tropospheric version accounts for a binary (H 2 O, 
H 2 SO 4 ) solution. 

The most prominent feature of the model is that it treats 
particle sizes in a Lagrangian scheme; i.e. the model follows 
the radii of growing/shrinking particles and creates a new 
size class i each time there is new ice or NAT particle nu- 
cleation. The model then transports the liquid, NAT and ice 
particles size class i downstream with time-dependent radius 
r,(t) and constant number density, 11 , . The Lagrangian size 
treatment avoids numerical diffusion, which may lead to ar- 
tificial redistributions of particle number densities on a fixed 
Eulerian size grid with fixed bins. 

The model is initialised with a log-normally distributed 
population of solution droplets, which represent the back- 
ground aerosol. This initial distribution is calculated given 
a user-defined mode radius, number density and standard de- 
viation, and in this study consisted of 26 radii. The initial 
composition of the droplets is calculated assuming thermo- 
dynamic equilibrium with the gas phase and accounting for 
the Kelvin effect of the various species. The total mixing ra- 
tios of H 2 O, H 2 SO 4 and, in the PSC version, HNO 3 , are set 
in the input file and are assumed to remain constant for the 
whole simulation along a trajectory (i.e. neglecting mixing 
and particle sedimentation). 

The model is run along a trajectory consisting of a time 
series of temperature and pressure data points, and after the 
initialisation, the aerosol droplet sizes are allowed to evolve 
freely in radius space, and are no longer restricted to a log- 
normal form (see Meilinger et al., 1995). As the temperature 
decreases, growth of the aerosol droplets is calculated, fully 
kinetically for all species, whereby the Kelvin effect is ac- 
counted for in calculating the vapour pressures of the indi- 
vidual species over the droplet surfaces. Latent heat effects 
of condensation and evaporation (or freezing) are not ac- 
counted for in the PSC version of the model; however given 
the small mixing ratios of matter which condense or evapo- 
rate (a few ppmv for H 2 O and a few ppbv for HNO 3 ) in the 
polar stratospheres, this does not impart any significant un- 
certainties to the calculations. Accommodation coefficients 
for the uptake/evaporation of species from/to the gas phase 
were set to 1 . 0 , for liquid and for solid particles. 

A schematic representation of the modelled composition 
and phase changes of PSCs as a function of temperature is 
shown in Fig. 2. Pathways previously represented in ZOMM 
are marked with dashed arrows; those which have been added 
for this study are shown by solid arrows. At temperatures 
above about 200 K, supercooled aerosol droplets exist as 
a quasi-binary H 2 O/H 2 SO 4 solution. Upon further cooling, 
HNO 3 is taken up to form a supercooled ternary solution 
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H,SOj / H,0 

T [K] liquid sulfate aerosol 



Fig. 2. A schematic description of the different pathways included 
in ZOMM. The dashed arrows depict processes which were rep- 
resented in previous model versions; the solid arrows show those 
added to the model for this study. 7d ew denotes the dew point, i.e. 
the temperature at which water saturation is reached. The tempera- 
tures indicated for 7d ew , 7f rost and 7 nat are approximate, based on 
typical polar stratospheric conditions. Foreign nuclei are indicated 
by the black triangles. While here the focus is on the heterogeneous 
nucleation of NAT (thick black arrow), the heterogeneous nucle- 
ation of ice (thin solid black arrows) is discussed in a companion 
paper (Engel et al., 2013). The heterogeneous nucleation of NAT 
occurs in STS aerosol droplets which contain a foreign nucleus (tri- 
angle surrounded by circle on the upper right of the figure), after 
the nucleation of NAT and as the air parcel continues to cool, a 
supercooled binary solution (SBS) coating, consisting of H 2 O and 
H 2 SO 4 , is present on the NAT particle. 

(STS). In the absence of solid nuclei, ice may nucleate homo- 
geneously in the STS droplets when temperatures decrease a 
few degrees below the frost point, and NAT can then form 
on the surface of the ice particles. Alternatively, if foreign 
nuclei are present (as indicated by the black triangles), they 
may lead to ice nucleation at higher temperatures than in the 
homogeneous case (Engel et al., 2013), or, as shown here, to 
the formation of NAT in the absence of ice (thick black line 
in Fig. 2). 

2.2.1 Homogeneous ice nucleation 

The homogeneous nucleation rate coefficient for the forma- 
tion of water ice, is calculated based on the tempera- 

ture and the water activity in the droplet according to 

•Ace™ = 10 D cm -3 s -1 , (1) 

where 

D = -906.688 + 8502.28* - 26 924.4.V 2 + 29 179.6* 3 (2) 

and * = fl w ,liq— «w,ice is the water activity in a droplet (« w ,ii q ) 
minus the water activity in a solution which would be in equi- 
librium with ice (a w ,ice) at the particular temperature being 


considered. « WJCC is defined as the ratio of the H 2 O vapour 
pressures over ice and pure water, and its value is given by 
Zobrist et al. (201 1), based on Koop et al. (2000), as 

fl w,ice = 

£XP ( 8.31441x7 )' (3) 

where T is the absolute temperature in Kelvin. The expected 
number density of ice particles (nice) nucleated from mod- 
elled liquid droplets with number density nij q and volume 
uiiq can be calculated by 

t 

nice = «liq(l- exp j 7 i 1 “ m uii q dr^. (4) 

to 

Instead of forming a new Lagrangian ice particle size class 
for every droplet size class, in every time step, which would 
be computationally impractical, the number of nucleated ice 
particles in each droplet size class is accumulated until this 
number exceeds either lxl0 -5 cm -3 , or half the number 
density of droplets in that size class, or half of the total num- 
ber density of ice particles in the model, whichever of the 
three criteria is satisfied first. 

After ice is nucleated in a liquid size class, the number of 
frozen droplets calculated using Eq. (4) is transferred from 
the liquid size class to the new ice size class. All ice particles 
in a particular size class thus have the same radius, and this 
radius increases as the ice particles take up water in compe- 
tition with the remaining liquid droplets. Upon warming, the 
evaporation/sublimation of species from the solid particles is 
calculated, until eventually a concentrated solution droplet is 
released and the number of particles is transferred back from 
the ice size class to the liquid size class from where it origi- 
nated. 

2.2.2 NAT formation 

Until now, NAT formation in ZOMM was only permitted via 
deposition nucleation on clean (uncoated, dry) surfaces of ice 
particles. The method used was thoroughly described by Luo 
et al. (2003); a summary of this description is provided here. 

Ice particles in the atmosphere are generally coated in a 
layer of liquid solution; however there may be dry patches 
on the particle, where the ice interacts directly with the gas 
phase. In rapid cooling events, the supersaturations with re- 
spect to NAT in the gas phase can grow to much larger val- 
ues than in the liquid due to non-equilibrium effects, favour- 
ing deposition nucleation on these exposed clean ice surfaces 
(Luo et al., 2003). The area of the dry patches was deter- 
mined by assuming a minimum thickness of the liquid layer 
of 28.5 x 10 -8 cm (approximately 10 monolayers; Luo et al., 
2003), and calculating the resulting surface coverage of the 
actual volume of liquid available to form the layer on the 
solid particles. The NAT nucleation rate coefficient on the 
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dry area of the ice particles was calculated based on the clas- 
sical nucleation theory for the heterogeneous nucleation of 
water ice, using Eq. (5): 




het,ice 


NAT 


(T) = 


kT 

h 


exp 


x ft exp 


AFdjffCD ' 

kT 

A G(T)f het (a) 
kT 


(5) 


where T is the temperature in Kelvin, k and h the Boltzmann 
and Planck constants respectively, n is the number density 
ofHNC>3 molecules on the ice surface (ca. 2 x 10 l4 cm -2 , at 
high NAT supersaturations; see Abbatt, 1997), AEdiff is the 
temperature-dependent diffusion activation energy which has 
to be surmounted for 2-D diffusion of HNO3 to take place on 
the surface of a preexisting foreign nucleus (here the ice), 
and G is the Gibbs free energy for the formation of a critical 
NAT embryo in the absence of the foreign nucleus, i.e. 


A G(T) = 


1 67 r 


2 3 
m er° 


( 6 ) 


3 [kT p IiiISnat )] 2 ' 

In Eq. (6), m is the molecular mass of NAT, a is the surface 
energy of NAT, and Snat is the NAT saturation ratio. The 
presence of the foreign nucleus reduces AG by the compati- 
bility factor /h et , defined as 

/het(a) = ^ (2 + cosa) (1 - cosa) 2 , (7) 

where a is the contact angle, a parameterisation of the hetero- 
geneous nucleation activity of the foreign nucleus (Seinfeld 
and Pandis, 2006). The possible values for a range between 
0°and 180°, leading to values of /h e t between 0 and 1; i.e. 
with a contact angle of 180° the foreign nucleus does not 
cause any enhancement in the nucleation rate. 

Now, as AFdiff is unknown, as are the compatibility factor 
/het and the surface energy of NAT, Luo et al. (2003) defined 
two new constants (denoted here as b and y), which they 
used to tune their nucleation rate coefficient equation, so that 
the model reproduced the observed NAT distribution: 

Tnat 06 = 6.24 x 10 24 cm -2 s -1 K -1 


T exp 


-b 


exp 


y (273.15) 3 


( 8 ) 


r 3 (ln(S NA T )) 2 

where the value 6.24 x 10 24 corresponds to kn/h. The values 
that Luo et al. (2003) suggested for these two constants were 


AFdjff 

b= — =2000K 

k 


(9) 


and 

AG(r)7’ 2 ln(5 N AT(T)) 2 /het(a) 3 

1 NAM j net v — 328 K 3 . (10) 

r 273.15 3 k ^ ; 


Note that in this parameterisation, it is implicitly assumed 
that the nucleation sites on the foreign nuclei are all of the 
same quality, as the compatibility factor is included in y and 
is therefore invariant. 


2.2.3 Extension to heterogeneous NAT nucleation and 
heterogeneous ice nucleation on foreign nuclei 


The microphysical box model has recently been extended so 
that in addition to homogeneous nucleation of water ice, and 
heterogeneous nucleation of NAT on the surface of ice par- 
ticles, it now includes the heterogeneous nucleation of ice 
on foreign nuclei as well as NAT formation on foreign nu- 
clei (both in the immersion mode). These two new nucle- 
ation pathways are based on the active site parameterisation 
of Marcolli et al. (2007). As foreign nuclei, a monodisperse 
population of solid particles with radii of 20 nm and with a 
number density of 7.5cm -3 is assumed. Measurements re- 
ported by Curtius et al. (2005) show that during January to 
March 2003 approximately 67% of aerosol particles in the 
Arctic polar vortex contained non-volatile residuals, while 
outside the vortex the value was much lower, at approxi- 
mately 24%. Similar values (between 30% and 70% de- 
pending on altitude above the tropopause) are reported for 
meteoritic particles in aerosol by Murphy et al. (2007), for 
mid- to low latitudes. During the RECONCILE (Reconcil- 
iation of essential process parameters for an enhanced pre- 
dictability of Arctic stratospheric ozone loss and its climate 
interactions) campaign, values of up to 80 % were found (von 
Hobe et al., 2013). The value used here, of 7.5 cm -3 aerosol 
particles containing foreign nuclei, corresponding to 50 % of 
the total aerosol, was chosen as a conservative estimate. As 
discussed below, only a small fraction of these nuclei actu- 
ally participate in the NAT nucleation; therefore there would 
only be a significant effect on the results presented here if a 
far lower percentage were to be chosen. 

According to the active site theory, freezing cannot be ini- 
tiated equally well over the entire surface of the foreign nu- 
cleus; rather it occurs preferentially at certain active sites. As 
a proxy for the quality of these sites, the contact angle may 
be used, which essentially parameterises the reduction in the 
nucleation energy barrier associated with a particular active 
site. In applying Eq. (5) to NAT nucleation on nuclei other 
than clean ice surfaces, the variations in contact angle for ac- 
tive sites of different quality must be accounted for. We there- 
fore define y', which differs from y of Luo et al. (2003) only 
in that it no longer includes the compatibility factor. The new 
parameterisation for the nucleation rate coefficient of NAT, 
on foreign nuclei immersed in a ternary H2O/HNO3/H2SO4 
solution, is therefore 


' / NAT° rei8n ( r ) = 6 - 24 x 10 24 cm -2 s -1 (T /A) exp 


—b 

~T~ 


y , (273.15) 3 / he t(a) 

r 3 (ln(S NAT )) 2 


( 11 ) 


where 5/ at is the saturation ratio of the STS with respect to 
NAT. 

In order to adjust the nucleation rate parameterisation so 
that modelled NAT distributions correspond to observations, 
the value of y' and the distribution of contact angles a for 
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the available active sites on the foreign nuclei must be deter- 
mined. Assuming that only the best active site on a particle 
is responsible for NAT nucleation, the occurrence probability 
of an active site with a particular contact angle on a particle 
can be described as 


Has (a) 


/— 51.0° 

Ppre X exp 

\a-«o 


) 


( 12 ) 


similarly to Marcolli et al. (2007). The parameter ao rep- 
resents the lowest contact angle considered, while Ppre is a 
constant (with units of deg -1 ). In the model, 1° contact an- 
gle bins were used, for contact angles starting from ao + 1°. 
The number density of foreign nuclei, containing at least one 
active site with a minimum contact angle a (i.e. the number 
of foreign nuclei in one particular contact angle bin), is then 
given by 


O'— 1 


— ( ^foreign ^ ' n a' 
\ a'=t*o+l 


Pas (a) Aa 


^foreign 


(13) 


where Aa = 1°, rf ore ign is the radius of the foreign nuclei 
(20 nm) and A\ is the area of one active site, which is taken 
to be 10 nm 2 , following Marcolli et al. (2007). The nucleation 
rate of NAT on a particle is then T^j forelgn multiplied by the 
area of the active site, A\. This corresponds to a volume- 
based nucleation rate of 


,het,vol, foreign 
•'NAT 


( T) = 


Q-max 

J n cli xJ^ oreiefl (T,a)xA 1 . (14) 
a,-=ao+l 


The number density of foreign nuclei as a function of the 
contact angle of their best active site is shown in Fig. 3. This 
figure highlights the contrast between the very few excel- 
lent nuclei, characterised by small contact angle active sites, 
and the large number of mediocre nuclei with larger con- 
tact angles. Similar distributions of nuclei quality have pre- 
viously been observed for populations of volcanic ash par- 
ticles (Hoyle et al., 2011) and mineral dust particles (Pinti 
et al., 2012) in ice nucleation experiments. The curve in 
Fig. 3 becomes rather flat beyond 80°, indicating there are 
similar numbers of particles having a best active site with 
a contact angle of, for example, 100° or 120°. With a for- 
eign nuclei number density of 7.5 cm -3 , and setting Pprc — 
1 x 10 -6 deg -1 (see below), considering contact angles up to 
146° accounts for approximately 0.13 particles cm -3 , which 
covers the range of NAT number densities necessary for sim- 
ulating the PSCs observed during the winter of 2009/2010. 
Note that Eq. (12) is not normalised, which effectively means 
that the majority of the foreign nuclei have best active sites 
which are of insufficient quality to decrease the NAT nucle- 
ation barrier. 

The properties of the PSC modelled with Eq. (11) can be 
adjusted by varying the values of the three fitting parameters 
y ' , ao and P pre . A higher value of y' effectively increases the 
nucleation barrier represented by the last term of Eq. (11), 



Contact angle [deg.] 


Fig. 3. The number density of foreign nuclei as a function of the 
contact angle of the best active site on the particle’s surface. 


and therefore lowers the temperature at which NAT will first 
form. The parameter ao sets a lower limit on the contact an- 
gle associated with the best possible active site. As described 
above, an active site with a small contact angle has a low 
value of the compatibility factor (/het); therefore ao deter- 
mines the maximum extent to which the presence of a foreign 
nucleus can reduce the nucleation barrier. On the other hand, 
Hpre changes the probability of an active site with a particu- 
lar contact angle occurring, and therefore alters the number 
density of NAT that can be formed at a particular tempera- 
ture. As it is assumed that the best active site on a foreign 
nucleus is responsible for the nucleation event, an increase in 
Ppre shifts the distribution of active sites towards lower con- 
tact angles, increasing the number of NAT particles which 
nucleate at higher temperatures. 

While for ice the nucleation rate is the sum of the hetero- 
geneous and homogeneous rates, in the case of NAT forma- 
tion only the heterogeneous nucleation rate on foreign nuclei 
and on ice is considered. Upon evaporation of the ice or NAT 
particles, the nuclei are returned to the active site bin from 
which they originated. 

The choice of the fitting variables for NAT nucleation, as 
well as the other parameters used in this particular study, is 
discussed in Sect. 2.3, while the parameters used for hetero- 
geneous ice nucleation are given in Engel et al. (2013). 

The optical properties (backscatter ratio, P 532 , and aerosol 
depolarisation ratio, <5 ae rosol) of the modelled PSC are cal- 
culated from the particle number densities and radii, as- 
suming the particles to be prolate spheroids with aspect ra- 
tios (diameter-to-length ratios) of 0.9 (in accordance with 
Carslaw et al., 1998) and a refractive index of 1.48 for NAT 
(Middlebrook et al., 1994). Calculations were then made us- 
ing a T matrix light scattering algorithm (Mishchenko, 1991; 
Carslaw et al., 1998) for spheroids with a known refractive 
index. 
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Fig. 4. Modelled aerosol depolarisation ratio (<5 a erosol) vs. inverse backscatter ratio ( l//?532) without (panel a) and with (panel b) consideration 
of instrument noise, for the mixed NAT/liquid cloud observed on the third CALIPSO overpass of the Arctic vortex on 21 December (denoted 
here as orbit 21_03). The error bars illustrate one standard deviation of the uncertainty for a single point. The different areas of the plots 
correspond to the PSC classifications described in Sect. 2.1. 


The CALIOP measurements of parallel and perpendicular 
backscatter are subject to an uncertainty, and to simulate this 
in the model results, a normally distributed random instru- 
mental noise component is added to each of these quantities 
in the model output, before calculating .R 532 and ^aerosol- Sub- 
sequently, the modelled PSC is classified according to the 
same scheme used for the observations. This allows a direct 
comparison between the modelled and observed PSC proper- 
ties. The calculation of the uncertainty in the CALIOP mea- 
surements is described in detail in Engel et al. (2013), and 
an example of the model output without and with the instru- 
mental noise is shown in Fig. 4, panels a and b respectively. 
Additionally, error bars show the limits of one standard devi- 
ation of the normally distributed noise component for a sin- 
gle modelled point. The main effect of the random noise is to 
increase the variability in the depolarisation ratio. 

2.3 Trajectory calculations 

The microphysical model was run along trajectories which 
ended at CALIPSO observation points. For every CALIPSO 
orbit during December 2009, where PSCs were observed, 
the following criteria were applied to choose the observation 
points used in the analysis: 

- Latitude north of 50° N; 

- Altitude ofbetween 18km and 26km; 

- The temperature at the observation point was 7 nat or 
less; 

- Every 5th vertical profile available in the CALIPSO 
data set was chosen, to give a 25 km horizontal resolu- 
tion; 


- Every vertical level was used, providing approximately 
a 180m vertical resolution. 

From all the points fulfilling these criteria, ten-day back- 
ward trajectories were calculated, using the trajectory mod- 
ule of CLaMS (McKenna et al., 2002). The trajectories were 
based on ECMWF ERA-Interim data (Dee et al., 2011) with 
a resolution of 1° x 1° and 6 h, using potential temperature as 
the vertical coordinate. The vertical resolution was approx- 
imately 1.25km between the altitudes of 19km and 30km 
at polar latitudes. The internal and output time step for the 
trajectory calculation was 15 min. To determine the starting 
point for the model calculation, the trajectories were fol- 
lowed backward in time from their end points, until the tem- 
perature rose above 7nat, so as to avoid the complication 
of having to initialise the PSCs. If no such point was en- 
countered within the ten days preceding the CALIPSO mea- 
surement on which the trajectory ended, the trajectory was 
not used in the modelling. For each of the valid trajectories, 
the model was initialised and run forward in time up until 
the point of the CALIPSO observation. For the days which 
were investigated in December, there were around 2000- 
5000 valid trajectories per orbit. 

In previous work on cirrus nucleation (e.g. Karcher and 
Lohmann, 2002; Hoyle et al., 2005), as well as in the anal- 
ysis of ice PSCs (Engel et al., 2013), it has been shown that 
the observed microphysical properties of a cloud can only 
be properly reproduced by a model if small-scale tempera- 
ture fluctuations (for example resulting from gravity waves 
with horizontal wavelengths of less than 100 km, i.e. below 
the resolution of CLaMS and the ERA-lnterim data) are ac- 
counted for. Such fluctuations are expected to be ubiquitous 
in the atmosphere (e.g. Gary, 2006, 2008) and have been the 
subject of numerous PSC studies in the past, most of which 
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examine the formation of ice and subsequent NAT nucleation 
on ice particles (e.g. Fueglistaler et al., 2002; Luo et al., 2003; 
Mann et al., 2005; Hopfner et al., 2006; Eckemiann et al., 
2009; Alexander et al., 2011). 

It was demonstrated by Meilinger et al. (1995) that, as the 
equilibration times for HNO 3 and for H 2 O in droplets differ 
greatly, during the warming phase of a temperature fluctua- 
tion, H 2 O may partition from a liquid droplet back into the 
gas phase much more quickly than HNO 3 . This can result 
in an enhanced NAT saturation ratio, and NAT nucleation 
at higher temperatures than would otherwise be possible in 
the absence of temperature fluctuations. Several sensitivity 
runs have been performed to test the effect of small-scale 
temperature fluctuations on NAT nucleation. The time series 
of temperature fluctuations which was superimposed on the 
CLaMS trajectories for this purpose was based on that de- 
rived in Hoyle et al. (2005), and had mean amplitudes of 
±0.5 K, similarly to those used by Engel et al. (2013), for 
heterogeneous ice nucleation. The results of these model runs 
are presented in Sect. 3.4. 

Model configuration used in this study 

For this study, the total number density of liquid aerosol par- 
ticles was set to 15 cm -3 . Half of these particles contained 
no foreign nuclei, and the initial radius of these particles fol- 
lowed a log-normal distribution, with mode radius of 70 nm, 
and a = 1.8 (e.g. Dye et al., 1992). The number density of 
foreign nuclei was set to 7.5 cm -3 , with a radius of 20nm. 
The foreign nuclei were immersed in the remaining 7.5 cm -3 
liquid droplets, which were assigned a radius of 70 nm (equal 
to the mode radius of the liquid-only droplet distribution). 
This was done to achieve the computational efficiency neces- 
sary to run the large number of trajectories simulated. As all 
NAT nucleation in the simulations occurred heterogeneously, 
this simplification had no significant effect on the results. 
The composition of the liquid solution was initially in equi- 
librium with the gas phase, and the H 2 O and HNO 3 mix- 
ing ratios were set at the beginning of each trajectory, us- 
ing daily mean values of Microwave Limb Sounder (MLS) 
measurements. These daily vortex average values were cal- 
culated from the MLS data by averaging over all cloud-free 
points within the polar vortex on a particular day. By adjust- 
ing the model until the CALIOP PSC observations were re- 
produced as well as possible, the following best values were 
determined for the fitting parameters: a, — 43°, y' = 650 K 3 
and P pre = 1 x 10 - 6 deg -1 . In Sect. 3, comparisons of the 
observations and model calculations using the new parame- 
terisation with these values are presented, as well as further 
calculations which illustrate the sensitivity of the results to 
changes in the fitting parameters. 


2.4 Constant NAT nucleation rate 

Experiments using constant, air-volume -based nucleation 
rates were also performed to provide a comparison with 
the new physically based NAT nucleation parameterisation. 
Similarly to previous studies (as summarised in Table 1), 
the constant nucleation rate was applied at all the trajec- 
tory points where the temperature was below 7nat- Rate co- 
efficients of between /^at ™ 1 = 7.0 x 10 -6 cm -3 airh -1 and 
^nat ™ 1 = ^ x 10 -5 cm -3 airh -1 were run along trajectories 
for the fourth CALIPSO overpass of the Arctic vortex on 24 
December 2009 (hereafter 24 04), and the nucleation rate 
producing the best results, /nat™* = 9 x 10 -6 cm -3 airh -1 , 
was then applied to trajectories for the remaining orbits. The 
results of these model experiments are also shown in Sect. 3. 


3 Results 

3.1 Constant vs. saturation-dependent nucleation rate 
coefficients 

The basic difference between the new parameterisation and 
the constant, volume-based nucleation rate can be illustrated 
by plotting the development of the NAT saturation ratio and 
the increase in NAT number density as a function of time, for 
a typical trajectory (taken from orbit 24 04) as in Fig. 5. The 
model was run once with each nucleation parameterisation, 
along this 60 h trajectory, which was initialised with a tem- 
perature of 193. 05K and a pressure of 25.5hPa. The HNO 3 
mixing ratio of 13.5ppbv and water vapour mixing ratio of 
5.7ppmv were taken from the daily average MLS data, as 
described in Sect. 2.3. 

The constant nucleation rate was applied as soon as the 
temperature drops below 7nat, and therefore the number 
density of NAT particles began to increase at this point. The 
new parameterisation on the other hand led to a significant 
formation of NAT only after the temperature reached approx- 
imately 7 nat — 4K. Shortly after this point, the tempera- 
ture stabilised briefly; therefore only about 1 x 10 -4 cm -3 of 
NAT was nucleated. Subsequently, the temperature contin- 
ued to decrease, and the number density of NAT nucleated 
with the new parameterisation climbed rapidly to approx- 
imately 1 x 10 -3 cm -3 . The step-like nucleation behaviour 
seen from the new parameterisation comes about because of 
the differing active site qualities on the particles. Once all 
the particles with sufficiently good active sites have nucle- 
ated NAT at a particular NAT saturation ratio, the NAT sat- 
uration ratio must increase before the particles with the next 
best active sites may also cause NAT nucleation. 

The saturation ratio in both simulations remained fairly 
similar while the air parcel was cooling; however due to the 
longer times available for the NAT particles to grow, after 
approximately 25 h the uptake of HNO 3 by the NAT parti- 
cles in the constant nucleation case led to a slightly lower 
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Fig. 5. A comparison of the volume-based nucleation rate parame- 
terisation (red) and the new parameterisation (black). Both versions 
of the model were run along the temperature trajectory shown in the 
bottom panel. The development over time, of NAT number density 
and saturation with respect to NAT, is shown in the top and middle 
panels respectively. 

saturation ratio than for the model with the new parameter- 
isation. After approximately 40 h, the NAT saturation ratio 
became lower for the simulation with the new parameterisa- 
tion, and continued to decrease faster than in the simulation 
with the constant nucleation rate. 

Three important differences between the constant nucle- 
ation rate and the new parameterisation are illustrated by 
Fig. 5. The first is that if the temperature is only slightly 
lower than 7 nat, the constant rate will nucleate NAT par- 
ticles, whereas the new parameterisation will not (0 to 20h 
in Fig. 5). Secondly, for short periods of time at low temper- 
atures, the new parameterisation will nucleate far more NAT 
than the constant rate (e.g. 20 to 3 0 h ) . Finally, if the temper- 
ature remains below Tnat for very long periods of time (as 
can be seen in Fig. 7, times below Tnat of 120h or more are 
not unusual), the constant rate will eventually nucleate more 
NAT than the new parameterisation. This would happen after 
approximately 1 10 h for the trajectory in Fig. 5, if it were to 
continue without further cooling. 

3.2 Sensitivity of new parameterisation to choice of 
parameters 

The values of the model parameters y' , ao and P prc were de- 
termined by systematically varying each parameter in turn, 
and comparing the resulting modelled PSC optical proper- 
ties with the CALIOP observations. For this purpose, the or- 
bit 24 04 was particularly suitable, as a fairly large, coherent 
area of PSC was observed, and the length of time that this 
air mass had experienced temperatures below Tnat was only 
between about 48 and 72 h, minimising complications caused 
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by mixing or by the sedimentation of PSC particles. After the 
match between the modelled results and the observations for 
orbit 24 04 was optimised, the model was applied to other 
orbits in a predictive capacity, to evaluate the robustness of 
the new parameterisation. 

The PSC classifications of the model runs used to define 
the model parameters are shown in Fig. 6, along with con- 
tours indicating the temperature at the observation point rel- 
ative to Tnat. The uppermost panel shows the observed PSC 
classifications. On the top row of the figure, the nucleation 
barrier parameter y' is varied from a value of 550 to 750 K 3 
in steps of 50 K 3 . As expected at low values of y', the PSC 
area is greater, and with increasing j / the area of PSC formed 
gradually decreases until only those trajectories that had ex- 
perienced the coldest temperatures still nucleate NAT. On the 
second row of Fig. 6, the nucleation strength parameter Tpre 
is varied from 1 x 10 -7 deg -1 to 1 x 10 -5 deg -1 . Here, the 
situation is reversed, with the larger values leading also to 
larger areas of PSC formed, and also an enhancement in the 
NAT number density, as shown by the increasing number of 
points being classed as Mixl/Mix2 and then as Mix-2enh. 
On the bottom row of the figure, the compatibility factor ao 
is varied from 41 to 45° in steps of 1° . Again, the increasing 
values of ao lead to a decrease in the temperature at which 
NAT can form, and a reduction of the area of NAT modelled. 
The three plots in the middle column of Fig. 6 are identical, 
and represent the model results using our chosen set of pa- 
rameters (a 0 = 43° , P pre = 1 x 10 -6 deg -1 and y' = 650 K 3 ). 

It is evident from this plot that the same results can be 
achieved with different combinations of the three factors. For 
example, an ao of 43°, P pre of 5 x 10 -7 deg -1 and y' of 
650 K 3 leads to a very similar PSC classification as is found 
when using an ao of 45°, P pre of 1 x 10 -6 deg -1 and y' of 
650 K 3 . We do not, therefore, claim that the values for the pa- 
rameters which we have determined are the only values that 
will result in an accurate modelling of the observed PSC. 
However, lacking laboratory data on quantities such as the 
distribution of the quality of the foreign nuclei (i.e. the active 
site distribution), A Pdiff and the surface energy of NAT, it 
is not possible to further constrain the choice of parameters, 
and the set of values which we have chosen are shown, in 
Figs. 7, 8 and 9, to lead to a very good representation of the 
observed NAT PSCs. Additionally, it can be seen in Fig. 6 
that the results are relatively insensitive to small changes in 
the three parameters, indicating a certain robustness of the 
parameterisation. 

In Fig. 7, PSC classifications calculated from observa- 
tional data (top row) and model output from simulations us- 
ing the new parameterisation (second row) as well as the con- 
stant, volume-based nucleation rate (third row) are shown. 
The plots represent a selection of orbits from throughout 
late December in which the advantage of the saturation- 
dependent nucleation rate parameterisation over a constant 
rate is particularly evident. The bottom two rows of Fig. 7 
show the time that the air parcel spent at temperatures below 
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Fig. 6. An illustration of the sensitivity of the modelled NAT classification (the different classifications are indicated by the colours in the 
plot) to the three fitting parameters used in the new model, applied to CALIPSO orbit 24-12-2009 05:22:08 UTC (denoted orbit 24_04 in 
the text). The uppermost panel shows the observed PSC, and the contour lines show the temperatures with respect to 7jyAT (from the ERA- 
Interim reanalysis) at the time of CALIOP observation. The top row of panels shows how the results change as y' is varied, while /p re is 
held constant at 1 x 10 -6 deg -1 , and ag at 43° . In the second row, /p re is varied while y' is held constant at 650K 3 , and <xq at 43° , while 
in the bottom row Pp re and y' are held constant at 1 x 10 -6 deg - 1 and 650 K 3 respectively, while t*o is varied. 


7nat and below Tnat - 4K, while the contour lines on these 
two rows of plots show the minimum temperature reached 
along the trajectory, with respect to That- 

Orbit 24_04. The middle column of Fig. 7 shows the re- 
sults for orbit 24 04, which, as described above, was used 
to fit the model parameterisations. The observed PSC con- 
sisted of a compact area of Mixl and Mix2, with a few spots 
of Mix2-enh and STS. The new parameterisation reproduces 
the observed area of PSC rather well; however Mix2 is over- 
represented, and there are only a few areas of Mixl about the 
periphery of the cloud. The constant nucleation rate produces 
considerably less PSCs than observed, and at a lower alti- 
tude than the observed PSC. Of course a larger area of PSC 
could have been formed with a higher constant rate; however 
other orbits shown in Fig. 7 illustrate that a higher rate will 
not improve the overall agreement between the constant rate 
simulations and the observed PSCs. 

Orbit 19_08. In the first column of Fig. 7, results for orbit 
19 08 are shown. These were among the first observations of 
PSC for the 2009/20 10 winter, and the observed area of Mixl 
and Mix2 PSC is small. The new parameterisation captures 
this well, also producing a small speck of PSC in the correct 


general area of the plot, while the constant nucleation rate, 
despite being on the low side for orbit 24_04, produces far 
too great an area of PSC. The reason for this can be seen by 
examining the temperature plots for this orbit. While a large 
number of trajectories spent up to approximately 100 h below 
?nat, only a very small fraction reached temperatures below 
?nat - 4 K. As seen in Fig. 5, while the constant nucleation 
rate begins to produce NAT at a slow rate as soon as the tem- 
perature is below Tnat, the new parameterisation requires 
higher saturations and, once these saturations are reached, 
produces NAT relatively quickly. Of the studies discussed 
above, which used an air-volume-based nucleation rate, only 
GrooB et al. (2005) carried out any simulations where a cer- 
tain NAT saturation ratio had to be exceeded before NAT was 
nucleated. In their case, a constant nucleation rate was ap- 
plied as soon as the temperature was below Tnat, to form 
a low number density mode of large, sedimenting particles, 
while a mode of higher number density and smaller, non- 
sedimenting particles was nucleated only once the NAT sat- 
uration ratio exceeded 30 (equivalent to approximately Tnat 
— 5K). Applying such a nucleation barrier may slightly im- 
prove the agreement of a simple volume-based nucleation 
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Fig. 7. A comparison of the PSC classifications calculated from CALIOP observations (top row), from model results using the new nucleation 
parameterisation (second row) and from model results using a constant, volume-based nucleation rate (third row). The contour lines in the 
upper three rows indicate the temperatures with respect to That (from the ERA-Interim reanalysis) at the time of CALIOP observation. Each 
column shows the results from one orbit. The bottom two rows of plots show the time spent by the trajectories at temperatures lower than 
T nat (fourth row) and TjvjAT — 4 K (bottom row). On these two rows, the contour lines show the minimum temperature encountered along 
the trajectory, with respect to 7 nat- 


rate with the observations in some situations; however it will 
not solve the fundamental problem that such a nucleation pa- 
rameterisation cannot reproduce both PSCs formed during 
short periods at cold temperatures and PSCs formed during 
long periods at relatively warm temperatures. 

Orbit 20_04. In the next column, results from the orbit 
20 04 are shown, where the new parameterisation this time 
does not do as well, underestimating the area of PSC formed. 
The results from the constant nucleation rate are, however, 
again an overestimate of the PSC area. Orbit 20 04 is rather 
unique, in that the time spent at low temperatures is very 
short compared to the other orbits. The observed PSCs are 
mostly STS and Mixl clouds, an aspect which is missing in 
the model results, which suggests that the ERA-Interim tem- 
peratures along this orbit are very likely too warm. 

Orbit 24_05. For this orbit, the new parameterisation pro- 
duces PSCs at the correct altitude and longitude in compari- 
son with the observed PSCs, albeit a slightly too small area, 
while the constant nucleation rate leads to a total absence of 


detectable PSC, due to the short time the trajectories spent at 
temperatures below Tnat- 

Orbit 26_04. For this orbit, the constant nucleation rate 
leads again to a far greater area of PSC than that which was 
observed. The new parameterisation simulates a slightly too 
large area of PSC, again with an overestimation of the frac- 
tion of Mix2 with respect to Mixl. While only a few pixels 
of PSCs were observed at higher longitudes, the new parame- 
terisation leads to a significant region of a Mixl/Mix2 cloud. 
However, these simulated trajectories have spent more than 
80 h below 7nat — 4 K before reaching the orbit path. As dis- 
cussed above, these particles will have sedimented and mixed 
into air parcels with different saturations ratios, which is not 
taken into account in the analysis. This area has therefore 
been marked by shading. 
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Fig. 8. PSC classification (upper row) and optical properties obtained by CALIOP along several orbits in December 2009. The contour 
lines in the upper row show the temperatures with respect to 7 nat (from the ERA-Interim reanalysis). The second row shows the aerosol 
backscatter ratio (R 532 — 1) and the third the perpendicular backscatter coefficient (/iperp)- In the panels on the bottom row, each observation 
point is represented by a dot, placed according to the observed aerosol depolarisation (^aerosol); an d inverse backscatter ratio, measured at a 
wavelength of 532nm ( 1 /^ 532 ). The points are colour coded to indicate the temperature with respect to 7 nat at the observation point. 


3.3 Modelling the NAT formation in December 2009 


The ability of the new NAT nucleation parameterisation to 
reproduce the observed PSC characteristics is further illus- 
trated by comparing the results from simulating a range of or- 
bits from mid- to late December 2009. In Fig. 8, each column 
shows the classification and observed optical characteristics 
of PSCs observed along a particular orbit. The correspond- 
ing model results are shown in Fig. 9. The model reproduces 
the PSC classes quite well in all orbits, both in terms of the 
location and the size of the modelled PSC area, and largely 
also in terms of optical characteristics. However, the model 
slightly underestimates the observed proportion of STS. This 
underestimation would be reduced if the temperatures were 
lower. 

In the plots of aerosol backscatter ratio (R532 — 1). the 
observations exhibit a larger variability than the modelled 
results, with the observational data containing some points 
reaching values of up to about 5-6 (i.e. I/R532 > 0.15), while 
the modelled values are generally below 1 (i.e. VR532 > 0.5). 


The observations show perpendicular backscatter coeffi- 
cients (/Iperp) up to 1 x 10 -5 km -1 sr -1 scattered throughout 
the observation area, while in the modelled case such high 
values are confined to the vicinity of PSCs, and the mod- 
elled background signal remains below 1 x 10 -6 km -1 sr -1 . 
The cause of this difference is that the instrumental noise 
spikes in the background observations, particularly those 
which form the tails of the /3 per p distribution, are not entirely 
reproduced by the function used to generate the instrumental 
noise which was added to the modelled optical values. The 
noise in the measured backscatter is often caused by vari- 
ations in radiation from space, to which the satellite-based 
instruments are subjected. 

Comparing the bottom panels of Figs. 8 and 9 shows that 
while the model generally reproduces the distribution of the 
points as a function of <5 a erosol and inverse backscatter ratio, 
the distribution of the modelled points is more compact with 
respect to the inverse backscatter ratio than is the case for 
the observational data. While the points at small backscatter 
ratios (I/.K 532 > 0.5) are well represented, the model misses 
the observed large backscatter ratios. This concerns mainly 
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Fig. 9. As for Fig. 8, except that in this case the model results are presented. 


STS clouds. Since STS clouds, owing to their large surface 
area, respond to temperature changes without significant ki- 
netic delay, they act as thermometer; i.e. the only reasonable 
explanation is that the ERA-lnterim reanalysis temperatures 
are occasionally too high. In the critical temperature regime 
of maximum HNO 3 uptake by STS, a temperature reduction 
of just 1 K can cause STS clouds to grow by up to a fac- 
tor of 4 in volume (Carslaw et al., 1994), i.e. about a factor 
1.6 in radius, which translates into an increase in (R 532 — 1 ) 
by up to a factor 6 (Biele et al., 2001). A sensitivity test, in 
which the temperature along the entire length of the trajec- 
tories was reduced by 1 K, led to inverse backscatter ratios 
which spanned the range of the observations, and also to an 
increased area of the modelled PSC (not shown). Increasing 
the total number of aerosol droplets in the model which did 
not contain foreign nuclei from 7.5 to 15 cm -3 , on the other 
hand, increases the proportion of PSC classified as STS, but 
does not significantly increase the modelled backscatter ra- 
tios. The uncertainty of the ERA-lnterim stratospheric tem- 
peratures is of the order of IK during the 2009/2010 Arctic 
winter, as a detailed comparison with radiosondes launched 
from Ny Alesund and Sodankyla reveals (during the REC- 
ONCILE and LAPBIAT-II (Lapland Atmosphere-Biosphere 
Facility) campaigns, not shown here). This corroborates the 


notion that the lack of high backscatters is most likely due to 
uncertainties in the reanalysis data. 

3.4 The effect of temperature fluctuations 

The model runs shown in Fig. 9 were repeated; however 
this time, random phases of the temperature fluctuations de- 
scribed in Sect. 2.3 were superimposed on the trajectories 
along which ZOMM was run. As expected, due to the in- 
creased NAT saturation ratios resulting from the temperature 
fluctuations, the area of the modelled NAT PSCs increased, 
as did the number density of the NAT particles in the clouds. 
In order to restore agreement with the observed properties of 
the PSCs, the value of y' was increased from 650 to 700K 3 , 
increasing the NAT nucleation barrier. The results of these 
model runs are shown in Fig. 10. Comparing Fig. 10 with 
Fig. 8 , it can be seen that the areas and classifications of PSCs 
modelled when accounting for temperature fluctuations gen- 
erally agree very well with the observed PSC properties, al- 
though for orbits 2 103 and 22_03 the modelled PSC areas 
in Fig. 10 are slightly larger than observed. However, includ- 
ing the fluctuations also leads to a better representation of the 
observed R 532 — 1 and ft [K , v values, particularly in the case 
of orbits 20 08 and 23 05. Because of the effect of temper- 
ature fluctuations on the modelled PSC classifications, it is 
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Fig. 10. As for Fig. 9, except that in this case the model results were calculated along trajectories with small-scale temperature fluctuations 
included, and y' was set to 700 K 3 . 


Table 2. The parameters to be used in the NAT nucleation param- 
eterisation, for models which include a representation of the atmo- 
spheric dynamics leading to small-scale temperature fluctuations, 
and for models which do not include such a representation. 


Model dynamics 

/ (K 3 ) 

Ppre (deg ! ) 

“0 

No fluctuations 

650 

1 X 10“ 6 

43° 

Including fluctuations 

700 

1 X 10 -6 

43° 


useful to define two sets of model parameters, as shown in 
Table 2. The first set are for use in models where small-scale 
air motions and the associated temperature fluctuations are 
not explicitly accounted for. These parameters reproduce the 
effect of the temperature fluctuations on NAT formation, to 
a large extent. The second set of parameters is suitable for 
models where the temperature fluctuations are explictly rep- 
resented. 


4 Summary and conclusion 

The 2009/2010 Arctic winter was unique, in that during the 
second half of December the stratospheric temperatures were 
below Tnat, while having stayed at least 2K above the frost 
point and 5K above the temperature required for homoge- 
neous ice nucleation. Whereas there have been previous stud- 
ies which have suggested that nucleation on ice is not the 
only NAT formation pathway in the polar stratosphere, the 
CALIOP measurements presented by Pitts et al. (2011) and 
modelled here provide unprecedented evidence for persistent 
ice-free NAT formation on the synoptic scale. Because of the 
meteorological situation, the CALIOP observational data set, 
together with a microphysical model, could be used to per- 
form a detailed modelling of NAT formation in the absence 
of ice, enabling the derivation of a parameterisation of the 
heterogeneous nucleation process. A new NAT nucleation 
parameterisation was constructed, whereby NAT nucleates 
on foreign nuclei, such as meteoritic dust particles, which 
are immersed in stratospheric aerosol droplets. The parame- 
terisation uses active site theory to account for the range of 
qualities of the nucleation sites on the available nuclei. It ex- 
tends the parameterisation used by Luo et al. (2003) for NAT 
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formation via deposition nucleation on the uncoated surface 
of ice particles. 

The PSC model ZOMM, containing the new parameteri- 
sation, was applied to ten orbits between 1 9 December 2009, 
when the first pixels of PSC were observed by CALIPSO, 
and 26 December 2009. Later orbits were excluded as they 
were getting close to the date when ice was first observed and 
also because the time that a substantial fraction of the trajec- 
tories spent below Tnat was longer than 80 h, increasing the 
influence of factors other than nucleation on the reproducibil- 
ity of the observed PSC. 

It has been shown that the new NAT nucleation parameter- 
isation leads to an accurate modelling of the observed PSC 
throughout this period, with the extent of the observed NAT 
PSC being extremely well reproduced. There is one notable 
exception, orbit 20 04, where the area of PSC was under- 
estimated by the model. In this particular case, there was 
a relatively small number of trajectories which had spent a 
short time (less than 24 h) at low temperatures (Tnat — 4 K). 
Under such conditions, small differences between the actual 
meteorological situation and the ERA-lnterim data used by 
the model along the trajectories may lead to a larger apparent 
discrepancy between the modelled and observed PSC classi- 
fications. An examination of the modelled backscatter shows 
that it just marginally misses the CALIOP PSC classification 
criteria. 

The classification of the modelled PSCs corresponds very 
well with the classification of the observed PSCs; however 
the model has a slight tendency to over-represent Mix2, and 
underestimate Mixl and STS. Again this could be due to 
slightly too high temperatures in the trajectories. On the other 
hand, Engel et al. (2013) found that using the ECMWF op- 
erational analysis to calculate the trajectories, rather than the 
ERA-lnterim data, would have led to temperatures along the 
trajectories having a low bias. Further, it should be noted 
that although the classification Mix2 is slightly over repre- 
sented in the model results, the inverse backscatter values 
derived from both the modelled and observed PSCs straddle 
the boundary between Mixl and Mix2 and are actually very 
similar (e.g. bottom rows of Figs. 8 and 9). Therefore, small 
shifts in inverse backscatter lead to different modelled and 
observed classifications of the PSC despite very similar opti- 
cal properties. 

The model also reproduces the individual optical quanti- 
ties such as the backscatter ratio (R 532 ) and the perpendicular 
backscatter coefficient (/Sp erp ) very well, while the variability 
as well as peak values in these quantities are not quite as high 
as those observed. The most significant differences are in the 
fiperp outside of the areas where PSC has been detected. In 
these places the observed j9 per p exhibits a substantially larger 
variability than the modelled values. However, as no PSCs 
were detected in these areas, this has no impact on our re- 
sults, and is most likely due to instrumental noise. 

Small-scale temperature fluctuations, resulting for exam- 
ple from gravity waves, with horizontal wavelengths that are 


not resolved in the ERA-lnterim data, are expected to be 
present in the polar stratospheres, and indeed have been the 
topic of many previous PSC studies (e.g. Fueglistaler et al., 
2002; Luo et al., 2003; Mann et al., 2005; Hopfner et al., 
2006; Eckermann et al., 2009; Alexander et al., 2011). When 
such fluctuations are added to the trajectories along which 
ZOMM is run, higher NAT saturation ratios are reached, 
and thus, on average, NAT nucleates at higher temperatures. 
However, a similar effect can also be achieved to a large ex- 
tent with a lower NAT nucleation barrier (y') in the new pa- 
rameterisation. The standard value of 650 K 3 which was de- 
termined for y' can therefore be viewed, to a certain extent, 
as including the main effects of the temperature fluctuations. 
On the other hand, if the new parameterisation were to be 
implemented in a model where the small-scale temperature 
fluctuations are resolved, a value of y' =700 K 3 should be 
used instead. 

The new parameterisation has significant advantages over 
an approach which applies a single, saturation-independent, 
air-volume-based nucleation rate as soon as the temperature 
drops below 7nat- Such a method cannot accurately model 
both situations with temperatures remaining slightly below 
Tnat for an extended period of time as well as situations 
with temperatures substantially below Tnat prevailing over 
only short time periods. The new parameterisation copes well 
with both situations (see for example Fig. 7, orbits 24 05 and 
26 04), while also partially reproducing the observed PSCs 
in more extreme cases (e.g. orbit 20_04). 

The new parameterisation is computationally efficient and 
simple to implement in large-scale models. Beyond the ap- 
plication in our microphysical box model without sedimen- 
tation, we hope that the application of this parameterisation 
in three-dimensional simulations of denitrification and ren- 
itrification in the polar stratospheric vortex will provide the 
opportunity to follow the development of the NAT particles 
over extended time periods, from nucleation through their 
growth and subsequent sedimentation, all the way to final 
evaporation or tropospheric removal. 
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